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We present preliminary results of a search for the decay B^ —f a^(980)7r~ among 22.7 million 

T(4S) -^ BB pairs collected by the BABAR detector at PEP-11. Using improved background sup- 

^ ' pression techniques and optimal signal extraction for rare decay searches, an excess of events over 

expected background is observed at the level of 3.7 standard deviations. This corresponds to the 



Qh! branching fraction B{B^ — > ag'(ag' -^ r]TT^)'ir ) = (6.2ll]2'5 =t 1.1) x 10 ®, where the first error is 



statistical and the second is systematic. The 90% confidence level upper limit is 11.5 x 10 



Submitted to the 
2Qth International Symposium on Lepton and Photon Interactions at High Energies 

7/23/2001—7/28/2001, Rome, Italy 



Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 
Work supported in part by Department of Energy contract DE-AC03-76SF00515. 



The BABAR Collaboration, 

B. Aubert, D. Boutigny, J.-M. Gaillard, A. Hicheur, Y. Karyotakis, J. P. Lccs, P. Robbe, V. Tisserand 
Laboratoire de Physique des Particules, F-74941 Annecy-le-Vieux, France 

A. Palano 
Universita di Bari, Dipartimento di Fisica and INFN, 1-70126 Bari, Italy 

G. P. Chen, J. C. Chen, N. D. Qi, G. Rong, P. Wang, Y. S. Zhu 
Institute of High Energy Physics, Beijing 100039, China 

G. Eigcn, P. L. Rcincrtscn, B. Stugu 
University of Bergen, Inst, of Physics, N-5007 Bergen, Norway 

B. Abbott, G. S. Abrams, A. W. Borgland, A. B. Brcon, D. N. Brown, J. Button-Shafcr, R. N. Cahn, 

A. R. Clark, M. S. Gill, A. V. Gritsan, Y. Groysman, R. G. Jacobsen, R. W. Kadcl, J. Kadyk, L. T. Kcrth, 

S. Kluth, Yu. G. Kolomcnsky, J. F. Krai, C. LeClerc, M. E. Levi, T. Liu, G. Lynch, A. B. Meyer, 

M. Momayezi, P. J. Oddonc, A. Pcrazzo, M. Pripstein, N. A. Roc, A. Romosan, M. T. Ronan, 

V. G. Shelkov, A. V. Tclnov, W. A. Wenzel 

Lawrence Berkeley National Laboratory and University of California, Berkeley, CA 94720, USA 

P. G. Bright-Thomas, T. J. Harrison, C. M. Hawkes, D. J. Knowlcs, S. W. O'Ncale, R. C. Penny, 

A. T. Watson, N. K. Watson 

University of Birmingham, Birmingham, B15 2TT, United Kingdom 

T. Deppermann, K. Goetzen, H. Koch, J. Krug, M. Kunze, B. Lewandowski, K. Peters, H. Schmuecker, 

M. Steinke 

Ruhr Universitdt Bochum, Institut fiir Experimentalphysik 1, D-44^80 Bochum, Cermany 

J. C. Andress, N. R. Barlow, W. Bhimji, N. Chevaher, P. J. Clark, W. N. Cottingham, N. De Groot, 
N. Dycc, B. Foster, J. D. McFall, D. Wallom, F. F. Wilson 

University of Bristol, Bristol BS8 1 TL, United Kingdom 

K. Abe, C. Hearty, T. S. Mattison, J. A. McKenna, D. Thiesscn 
University of British Columbia, Vancouver, BC, Canada V6T IZl 

S. Jolly, A. K. McKemey, J. Tinslay 
Brunei University, Uxbridge, Middlesex UBS 3PH, United Kingdom 

V. E. Bhnov, A. D. Bukin, D. A. Bukin, A. R. Buzykaev, V. B. Golubev, V. N. Ivanchenko, A. A. Korol, 
E. A. Kravchenko, A. P. Onuchin, A. A. Salnikov, S. L Serednyakov, Yu. L Skovpen, V. L Telnov, 

A. N. Yushkov 

Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

D. Best, A. J. Lankford, M. Mandelkern, S. McMahon, D. P. Stoker 
University of California at Irvine, Irvine, CA 92697, USA 

A. Ahsan, K. Arisaka, C. Buchanan, S. Chun 
University of California at Los Angeles, Los Angeles, CA 90024, USA 



J. G. Branson, D. B. MacFarlane, S. Prell, Sh. Rahatlou, G. Raven, V. Sharma 
University of California at San Diego, La Jolla, CA 92093, USA 

C. Campagnari, B. Dahmcs, P. A. Hart, N. Kuznetsova, S. L. Levy, O. Long, A. Lu, J. D. Richman, 

W. Verkerke, M. Witherell, S. Yellin 

University of California at Santa Barbara, Santa Barbara, CA 93106, USA 

J. Beringer, D. E. Dorfan, A. M. Eisner, A. Prey, A. A. Grillo, M. Grothe, C. A. Heusch, R. P. Johnson, 

W. Kroeger, W. S. Lockman, T. PuUiam, H. Sadrozinski, T. Schalk, R. E. Schmitz, B. A. Schunim, 

A. Seiden, M. Turri, W. Walkowiak, D. C. Williams, M. G. Wilson 

University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, CA 95064, USA 

E. Chen, G. P. Dubois-Felsmann, A. Dvoretskii, D. G. HitHn, S. Metzler, J. Oyang, F. C. Porter, A. Ryd, 

A. Samuel, M. Weaver, S. Yang, R. Y. Zhu 

California Institute of Technology, Pasadena, CA 91125, USA 

S. Devmal, T. L. Geld, S. Jayatilleke, G. Mancinelh, B. T. Meadows, M. D. Sokoloff 
University of Cincinnati, Cincinnati, OH 45221, USA 

T. Barillari, P. Bloom, M. O. Dima, S. Fahey, W. T. Ford, D. R. Johnson, U. Nauenberg, A. Olivas, 
H. Park, P. Rankin, J. Roy, S. Sen, J. G. Smith, W. C. van Hoek, D. L. Wagner 

University of Colorado, Boulder, CO 80309, USA 

J. Blouw, J. L. Harton, M. Krishnamurthy, A. Soffer, W. H. Toki, R. J. Wilson, J. Zhang 
Colorado State University, Fort Collins, CO 80523, USA 

T. Brandt, J. Brose, T. Colberg, G. Dahlinger, M. Dickopp, R. S. Dubitzky, A. Hauke, E. Maly, 
R. Mliller-Pfefferkorn, S. Otto, K. R. Schubert, R. Schwierz, B. Spaan, L. Wilden 

Technische Universitdt Dresden, Institut fiir Kern- und Teilchenphysik, D-01062, Dresden, Cermany 

L. Behr, D. Bernard, G. R. Bonneaud, F. Brochard, J. Cohen- Tanugi, S. Ferrag, E. Roussot, S. T'Jampens, 

Ch. Thiebaux, G. Vasileiadis, M. Verderi 

Ecole Polytechnique, F- 9 1128 Palaiseau, France 

A. Anjomshoaa, R. Bcrnct, A. Khan, D. Lavin, F. Muheim, S. Playfer, J. E. Swain 
University of Edinburgh, Edinburgh EII9 3JZ, United Kingdom 

M. Falbo 
Elon University, Elon University, NC 27244-2010, USA 

C. Borean, C. Bozzi, S. Dittongo, M. Folegani, L. Piemontese 
Universitd di Ferrara, Dipartimento di Fisica and INFN, 1-44100 Ferrara, Italy 

E. TreadweU 
Florida A&M University, Tallahassee, FL 32307, USA 

F. Anulli,[| R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, D. Falciai, G. Finocchiaro, P. Patteri, 

L M. Peruzzi,^! M. Piccolo, Y. Xie, A. Zallo 

Laboratori Nazionali di Frascati dellTNFN, 1-00044 Frascati, Italy 



^ Also with Universita di Perugia, 1-06100 Perugia, Italy 



S. Bagnasco, A. Buzzo, R. Contri, G. Crosctti, P. Fabbricatore, S. Farinon, M. Lo Vetere, M. Macri, 
M. R. Monge, R. Musenich, M. Pallavicini, R. Parodi, S. Passaggio, F. C. Pastore, C. Patrignani, 

M. G. Pia, C. Priano, E. Robutti, A. Santroni 

Universita di Genova, Dipartimento di Fisica and INFN, I-I6I46 Genova, Italy 

M. Morii 
Harvard University, Cambridge, MA 02138, USA 

R. Bartoldus, T. Dignan, R. Hamilton, U. Mallik 
University of Iowa, Iowa City, lA 52242, USA 

J. Cochran, H. B. Crawley, P.-A. Fischer, J. Lamsa, W. T. Meyer, E. I. Rosenberg 
Iowa State University, Ames, lA 50011-3160, USA 

M. Benkebil, G. Grosdidier, C. Hast, A. Hocker, H. M. Lacker, S. Laplace, V. Lepeltier, A. M. Lutz, 
S. Plaszczynski, M. H. Schune, S. Trincaz-Duvoid, A. Valassi, G. Wormscr 

Laboratoire de lAccelerateur Lineaire, F-91898 Orsay, France 

R. M. Bionta, V. Brigljevic , D. J. Lange, M. Mugge, X. Shi, K. van Bibber, T. J. Wenaus, D. M. Wright, 

C. R. Wuest 

Lawrence Livermore National Laboratory, Livermore, CA 94550, USA 

M. Carroll, J. R. Fry, E. Gabathuler, R. Garnet, M. George, M. Kay, D. J. Payne, R. J. Sloane, 

C. Touramanis 

University of Liverpool, Liverpool L69 3BX, United Kingdom 

M. L. Aspinwall, D. A. Bowerman, P. D. Dauncey, U. Egedc, L Eschrich, N. J. W. Gunawardane, 

J. A. Nash, P. Sanders, D. Smith 

University of London, Imperial College, London, SW7 2BW, United Kingdom 

D. E. Azzopardi, J. J. Back, P. Dixon, P. F. Harrison, R. J. L. Potter, H. W. Shorthouse, P. Strother, 

P. B. Vidal, M. L Wilhams 

Queen Mary, University of London, El 4^S, United Kingdom 

G. Cowan, S. George, M. G. Green, A. Kurup, C. E. Marker, P. McGrath, T. R. McMahon, S. Ricciardi, 

F. Salvatore, L Scott, G. Vaitsas 

University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United 

Kingdom 

D. Brown, C. L. Davis 
University of Louisville, Louisville, KY 40292, USA 

J. AlHson, R. J. Barlow, J. T. Boyd, A. C. Forti, J. FuUwood, F. Jackson, G. D. LafFerty, N. Savvas, 

E. T. Simopoulos, J. H. Weatherall 

University of Manchester, Manchester M13 9PL, United Kingdom 

A. Farbin, A. Jawahery, V. Lillard, J. Olsen, D. A. Roberts, J. R. Schieck 
University of Maryland, College Park, MD 20742, USA 

G. Blaylock, C. DaUapiccola, K. T. Flood, S. S. Hertzbach, R. Kofler, T. B. Moore, H. Staengle, S. Willocq 

University of Massachusetts, Amherst, MA 01003, USA 



B. Brau, R. Cowan, G. ScioUa, F. Taylor, R. K. Yamamoto 
Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, MA 02139, USA 

M. Milek, P. M. Patel, J. Trischuk 
McGill University, Montreal, Canada QC H3A 2T8 

F. Lanni, F. Palombo 

Universitd di Milano, Dipartimento di Fisica and INFN, 1-20133 Milano, Italy 

J. M. Bauer, M. Booke, L. Cremaldi, V. Eschcnburg, R. Kroeger, J. Reidy, D. A. Sanders, D. J. Summers 

University of Mississippi, University, MS 38677, USA 

J. P. Martin, J. Y. Nief, R. Seitz, P. Taras, A. Woch, V. Zacek 
Universite de Montreal, Laboratoire Rene J. A. Levesque, Montreal, Canada QC II3C 3J7 

H. Nicholson, C. S. Sutton 
Mount Holyoke College, South Hadley, MA 01075, USA 

C. Cartaro, N. CavalloFI G. De Nardo, F. Fabozzi, C. Gatto, L. Lista, P. Paolucci, D. Piccolo, C. Sciacca 
Universitd di Napoli Federico II, Dipartimento di Scienze Fisiche and INFN, 1-80126, Napoli, Italy 

J. M. LoSecco 
University of Notre Dame, Notre Dame, IN 46556, USA 

J. R. G. AlsmiUer, T. A. Gabriel, T. Handler 
Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA 

J. Brau, R. Frey, M. Iwasaki, N. B. Sinev, D. Strom 
University of Oregon, Eugene, OR 97403, USA 

F. Colecchia, F. Dal Corso, A. Dorigo, F. Galeazzi, M. Margoni, G. Michelon, M. Morandin, M. Posocco, 
M. Rotondo, F. Simonetto, R. Stroili, E. Torassa, C. Voci 

Universitd di Padova, Dipartimento di Fisica and INFN, 1-35131 Padova, Italy 

M. Benayoun, H. Briand, J. Chauveau, P. David, Ch. de la Vaissiere, L. Del Buono, O. Hamon, F. Le 
Diberder, Ph. Leruste, J. Lory, L. Rods, J. Stark, S. Versille 

Universites Paris VI et VII, Lab de Physique Nucleaire H. E., F-75252 Paris, France 

P. F. Manfredi, V. Re, V. Speziali 
Universitd di Pavia, Dipartimento di Elettronica and INFN, 1-27100 Pavia, Italy 

E. D. Frank, L. Gladncy, Q. H. Guo, J. H. Panetta 
University of Pennsylvania, Philadelphia, PA 19104, USA 

C. Angelini, G. Batignani, S. Bettarini, M. Bondioli, M. CarpincUi, F. Forti, M. A. Giorgi, A. Lusiani, 
F. Martinez- Vidal, M. Morganti, N. Neri, E. Paoloni, M. Rama, G. Rizzo, F. Sandrelli, G. Simi, 

G. Triggiani, J. Walsh 

Universitd di Pisa, Scuola Normale Superiore and INFN, 1-56010 Pisa, Italy 



Also with Universita della Basilicata, 1-85100 Potenza, Italy 



M. Haire, D. Judd, K. Paick, L. Turnbull, D. E. Wagoner 
Prairie View A&M University, Prairie View, TX 7744^, USA 

J. Albert, C. Bula, P. Elmer, C. Lu, K. T. McDonald, V. Miftakov, S. F. Schaffner, A. J. S. Smith, 

A. Tumanov, E. W. Varnes 

Princeton University, Princeton, NJ 08544 ^ USA 

G. Cavoto, D. del Re, R. Faccini,^ F. Ferrarotto, F. Ferroni, K. Fratini, E. Lamanna, E. Leonardi, 
M. A. Mazzoni, S. Morganti, G. Piredda, F. Safai Tehran!, M. Serra, C. Voena 

Universita di Roma La Sapienza, Dipartimento di Fisica and INFN, 1-00185 Roma, Italy 

S. Christ, R. Waldi 
Universitdt Rostock, D- 18051 Rostock, Germany 

P. F. Jacques, M. Kalelkar, R. J. Piano 
Rutgers University, New Brunswick, NJ 08903, USA 

T. Adye, B. Franek, N. I. Gcddcs, G. P. Gopal, S. M. Xella 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

R. Aleksan, G. De Domenico, S. Emery, A. Gaidot, S. F. Ganzhur, P.-F. Giraud, G. Hamel de 
Monchenault, W. Kozanecki, M. Langer, G. W. London, B. Mayer, B. Serfass, G. Vasseur, Ch. Yeche, 

M. Zito 

DAPNIA, Commissariat a I'Energie Atomique/Saclay, F-91191 Cif-sur-Yvette, France 

N. Copty, M. V. Purohit, H. Singh, F. X. Yumiceva 
University of South Carolina, Columbia, SC 29208, USA 

I. Adam, P. L. Anthony, D. Aston, K. Baird, J. P. Berger, E. Bloom, A. M. Boyarski, F. Bulos, 

G. Calderini, R. Glaus, M. R. Convery, D. P. Coupal, D. H. Coward, J. Dorfan, M. Doser, W. Dunwoodie, 

R. C. Field, T. Glanzman, G. L. Godfrey, S. J. Gowdy, P. Grosso, T. Himel, T. Hryn'ova, M. E. Buffer, 

W. R. Innes, C. P. Jessop, M. H. Kelsey, P. Kim, M. L. Kocian, U. Langenegger, D. W. G. S. Leith, 

S. Luitz, V. Luth, H. L. Lynch, H. Marsiske, S. Menkc, R. Messner, K. C. Moffeit, R. Mount, D. R. MuUcr, 

C. P. O'Grady, M. Perl, S. Petrak, H. Quinn, B. N. Ratcliff, S. H. Robertson, L. S. Rochester, 

A. Roodman, T. Schietinger, R. H. Schindler, J. Schwicning, V. V. Serbo, A. Snyder, A. Soha, 

S. M. Spanier, J. Stelzer, D. Su, M. K. Sullivan, H. A. Tanaka, J. Va'vra, S. R. Wagner, 

A. J. R. Weinstein, W. J. Wisniewski, D. H. Wright, C. C. Young 

Stanford Linear Accelerator Center, Stanford, CA 94309, USA 

P. R. Burchat, C. H. Cheng, D. Kirkby, T. I. Meyer, C. Roat 
Stanford University, Stanford, CA 94305-4060, USA 

R. Henderson 
TRIUMF, Vancouver, BC, Canada V6T 2A3 

W. Bugg, H. Cohn, A. W. Weidemann 
University of Tennessee, Knoxville, TN 37996, USA 



Also with University of California at San Diego, La JoUa, CA 92093, USA 



J. M. Izcn, I. Kitayama, X. C. Lou, M. Turcotte 
University of Texas at Dallas, Richardson, TX 75083, USA 

F. Bianchi, M. Bona, B. Di Girolamo, D. Gamba, A. Smol, D. Zanin 
Universita di Torino, Dipartimento di Fisica Sperimentale and INFN, 1-10125 Torino, Italy 

L. Bosisio, G. Delia Ricca, L. Lanceri, A. Pompili, P. Poropat, M. Prest, E. Vallazza, G. Vuagnin 
Universita di Trieste, Dipartimento di Fisica and INFN, 1-34127 Trieste, Italy 

R. S. Panvini 

Vanderhilt University, Nashville, TN 37235, USA 

C. M. Brown, A. Dc Silva, R. Kowalcwski, J. M. Roncy 
University of Victoria, Victoria, BC, Canada V8W 3P6 

H. R. Band, E. Charles, S. Dasu, F. Di Lodovico, A. M. Eichcnbaum, H. Hu, J. R. Johnson, R. Liu, 
J. Nielsen, Y. Pan, R. Prepost, L J. Scott, S. J. Sekula, J. H. von Wimmcrsperg-Toellcr, S. L. Wu, Z. Yu, 

H. Zobernig 

University of Wisconsin, Madison, WI 53706, USA 

T. M. B. Kordich, H. Neal 
Yale University, New Haven, CT 06511, USA 



1 Introduction 

We present a search for the decay B^ -^ aQ'(980)7r~, where the oq resonance^ is observed in its 
dominant decay channel aj" — > r/vr"''. The interest in this mode stems from its potential use for 
measuring the CKM angle a [j^, ^. It was pointed out in Ref. Q that, within the factorization 
assumption, the main tree contributions to the aovr decay amplitude vanish: they would imply 
forbidden second class currents. This simplifies the two-body analysis for the extraction of a, and 
can lead to an enhanced direct CP violation. 

After a preselection, we refine the background suppression using Multivariate Analyzer (MVA) 
tools (Neural Net (NN) and Fisher discriminants). We use both cut-based and shape-based analy- 
ses, the latter employing a maximum likelihood technique, which together provide complementary 
results and cross-checks. The cut optimization procedure used in the analysis does not rely on 
a prior branching fraction estimate. Both inclusive and exclusive control samples are used for 
systematics checks, particularly with regard to the rj and oq resonances P). 

2 The BaBAR Detector and the Dataset 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring, 
located at the Stanford Linear Accelerator Center (SLAC). PEP-II is an asymmetric e^e~ collider 
with a center-of-mass energy equal to the T(4S') mass. From Nov. 1999 to Oct. 2000, a total 
of 22.7 million BB pairs has been collected by BABAR, corresponding to an integrated on-peak 
luminosity of approximately 20.7 fb~ . In addition, 2.6 fb^ of off-peak data were taken during the 
same period: they have been used to validate the contribution to backgrounds resulting from e+e" 
annihilation into light qq pairs. 

The BABAR detector and its performance are described in Ref. Q. The innermost component 
consists of a 5-layer Silicon Vertex Tracker (SVT), providing the positions of charged tracks in the 
neighbourhood of the beam interaction point. It is followed by a 40-layer central Drift Chamber 
(DCH), immersed in a 1.5-T magnetic field, measuring the track momenta and providing a mea- 
surement of the specific ionization loss (dB/dx) for particle identification (PID). The main PID 
device is a unique, internally refiecting ring imaging Cherenkov detector (DIRC), covering the cen- 
tral region of BABAR. A Cherenkov angle K/tt separation of better than 4 standard deviations is 
achieved for tracks below 3 GeV/c momentum. Photons are detected by a CsI(Tl) electromagnetic 
calorimeter (EMC), which provides excellent angular and energy resolution with high efficiency 
for energy deposits above 20MeV. A superconducting solenoid, located around the EMC is itself 
surrounded by an iron fiux return, instrumented with Resistive Plate Chambers to identify muons. 

3 Analysis Method 

3.1 Candidate Selection 

A B^ —^ flg vr^ candidate contains a pair of oppositely-charged pions and two photons. Charged 
tracks are required to satisfy a set of track quality criteria, which includes cuts on their momenta 
(less than lOGeV/c), transverse momenta (above lOOMeV/c), and on the number of DCH hits (at 
least 20). The tracks are also required to originate in the vicinity of the beam-beam interaction 
point. Pion candidates must fail electron selection criteria and are required to have a momentum 

^ Charge conjugation is implied throughout this document, and the a^(980) resonance is denoted oq. 



in the center-of-mass (CM) above 2 GeV/c. A B candidate is rejected if the track not used to 
form the uq has a DIRC Cherenkov angle consistent with a kaon. Photons are identified as energy 
deposits in the EMC, unassociated with charged tracks. They are required to have an energy above 
80MeV in the laboratory frame (LAB), and must satisfy photon-like shower profile criteria. To 
be associated with an rj decay, a pair of candidate photons is required to have an invariant mass 
0.470 < m(77) < 0.615 GeV/c^, and the rj CM momentum must be larger than 0.9 GeV/c. The pion 
track and r/ candidate form an oq candidate if their invariant mass falls in the range 0.9 < m{r]TT) < 
1.08GeV/c^. Reconstruction of B^ candidates is done by vertexing all combinations of vr+Tr^ry 
candidates in each event and applying a quality requirement on the vr+vr" vertex. A B^ candidate is 



characterized by two kinematic variables: the beam energy-constrained mass rriEC = V -^beam ~ ?%' 
where -Ebeam is half the CM energy and pB is obtained by applying kinematic constraints to the 
four-momenta of the B daughters; and \AE\ = \Eb — -Ebcaml with Eb being the CM energy of the 
B candidate. A candidate is retained if rriEc > 5.21GeV/c^ and \AE\ < 0.25 GeV. 

3.2 Background Suppression 

Charmless hadronic modes suffer from large amounts of background from random combinations of 
tracks, mostly from light quark production. In the CM frame, this background typically exhibits 
a two-jet structure in contrast to the spherically symmetric T{4S) — > BB events. Efficient back- 
ground rejection is obtained by requiring the angle 9t between the thrust axis of the B candidate 
and the thrust axis of the rest of the event (ROE) to satisfy | cos^tI < 0.9. Denoting | cos^tpI as 
the minimum cosine of the two angles formed by the two most energetic tracks (or neutral clusters) 
with respect to the thrust axis of the event, we require | cos^tpI < 0.88. Combinatorial background 
within a candidate event arises mainly from low-energy photons. Compared with decay modes con- 
taining vr^'s, this is a minor concern for uott due to the higher r] mass: the fraction of events with 
more than one photon pair combination passing the selection cuts is at the percent level. For events 
with multiple candidates the one with the most energetic low-energy photon is retained. 

Further discrimination is achieved by combining event shape variables using MVA techniques. 
A common approach uses the Fisher discriminant (denoted standard Fisher below) proposed by the 
CLEO Collaboration |0] . This standard Fisher is built as a linear combination of the cosine of the 
angle between the B candidate momentum and the beam axis, the cosine of the angle between the 
B candidate thrust axis and the beam axis, together with nine energy deposits, each defined to be 
the energy of charged tracks and neutral clusters of the ROE whose directions are contained in nine 
concentric cones centered around the B direction. In our analysis, the choice of variables has been 
reconsidered taking into account the separation powerg, the correlations between the variables and 
the signal efficiency at fixed background rejection. For this purpose, 12 variables are selected and 
combined using either a linear (Fisher) or a non-linear (NN) MVA. Table || lists the variables which 
are retained together with their respective Fisher coefficients. Among thse are the six variables 



'^ The separation power of the normahzed signal and background distributions, S{x) and B(x), of a discriminating 
variable x, is defined by 



2 J S{x)+B{x) 




-2-10 1 2 

Fisher (new) output 



0.5 -0.25 0.25 0.5 
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Figure 1: Fisher (left-hand plots) and Neural Net discriminants (right-hand plots). The upper 
plots show the distributions for signal and background, where the signal is taken from Monte Carlo 
simulation and the background from on-peak sideband data. The points with error bars show the 
off-peak data, containing no signal. The lower plots show the signal and background efficiencies as 
a function of the cut applied (the small dots indicate the off-peak data). 



Lf\ L^^ (j = 0,2,6) defined by 



] ' J 



^(c,n) 



J.I-? 



2J Pi X I COS (/J I , 

«(c,n) 



(1) 



which are the momentum-weighted sums of the cosines of the angles between the ROE charged 
tracks {Lj) or neutral clusters (Lj) and the thrust axis of the B candidate. These variables 
provide a generalization of the discrete cones used in the standard Fisher discriminant. 

Figure ffl shows the distributions of the Fisher and NN discriminants for signal and background, 
the former taken from Monte Carlo simulation and the latter from on-peak sideband data (his- 
tograms) and from off-peak data (points with error bars). The lower plots show the resulting signal 
and background efficiencies as a function of a cut applied on the output values of the discriminants. 
Figure |2| depicts the background versus signal efficiencies for the standard Fisher discriminant and 
the 12-variable MVAs adopted in this analysis. A 17% relative increase of the signal efficiency is 
obtained for NN with respect to the standard Fisher at the benchmark of 5% background retention. 
Table g summarizes the performances of the three MVA discriminants. The NN, being the dis- 
criminant which provides the best signal efficiency, is used for the signal extraction in this analysis, 
while keeping the Fisher for cross-checks. 
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Figure 2: Comparison of the background versus signal efficiencies for NN (solid line), the Fisher 
(dashed line), and the standard Fisher (dashed-dotted line) discriminants. 



3.3 Cut-Based Analysis 

The signal extraction is based on counting data events found in a signal region (SR) and subtracting 
from this the expected background yield. To find an optimal set of signal selection criteria, balancing 
between background rejection and signal efficiency, we use the procedure described in the appendix. 
The criterion it applies is to minimize the expected confidence level for the null hypothesis "there 
is no signal in the data sample" . Cut optimization is performed on the MVA output only, resulting 
in the requirements x°pj;(NN) > 0.6 and 2;°"]: (Fisher) > 0.7, respectively. 

The SR is defined as \m^c - 5.28| < 0.006 GeV/c^, |A£;| < 0.070 GeV and 0.517 < m(77) < 
0.587 GeV/c^. The background contamination in the SR is estimated from the two-dimensional 
Grand Sideband (GSB): |A^| < 0.25 GeV and 5.210 < m^c < 5.263 GeV/ c^, assuming for the 
corresponding background shapes a second-order polynomial for A£^, and an ARGUS function ||8[ 
for mEC- The signal efficiencies (e), event yields {M) and the expected backgrounds (A^^) are 
given in Table ^ for Fisher and NN. The corresponding branching fractions are obtained from the 
expression 



B 



eN 



BB 



(2) 



where N^j^ = (22.74 ± 0.36) x 10^ is the total number of produced BB pairs. Equation || assumes 
an equal production of charged and neutral S's in T{4S) decays. Non-resonant contributions to 
the r]7r~^Tr~ final state are disregarded for all branching fractions quoted in this document. The 
results obtained using Fisher and NN are compatible. Their statistical significances, defined as 
the probability to observe the given excess when no signal is present, are 3.0 and 3.1 standard 
deviations, respectively. The systematic errors assigned to the branching fractions are described in 
Sec. H. Also given in Table H are the upper limits at 90% confidence levels, where the systematic 
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Table 1: The 12 event-shape variables used as MVA inputs with their corresponding Fisher 
coefficients. 



Variable name Description 



Fisher coefficient 



R2 

I cos 0T I 
I COS ^S I 

I cos{B,z)\ 

|cos(B(T),z) 
I COS 6'tp I 

^2 

-^6 
re 



re 



Second Fox- Wolfram moment 
Angle: i?/ROE thrust axis 
Angle: B/ROE sphericity axis 
Angle: B direction/beam axis 
Angle: B thrust/beam axis 
Minimum cosine of ROE tracks/clusters 
Neutral zeroth-order angular function 
Neutral second-order angular function 
Neutral sixth-order angular function 
Charged zeroth-order angular function 
Charged second-order angular function 
Charged sixth order angular function 



0.40 
-0.71 
-0.89 
-0.88 
-0.79 

0.54 

0.35 
-1.02 
-0.69 

0.38 
-0.51 
-0.66 



Offset 



Centers the sum of signal and background at zero 



1.49 



Table 2: Results on the separation and the signal efficiency at the benchmark of 5% background 
retention for the standard Fisher discriminant and the two MVAs used in this analysis. 



Type 



::t\ 



e(5)(at e{B) = 5%) 



Fisher (standard) 0.44 

Fisher (12-var.) 0.48 

Neural Net (12-var.) 0.50 



0.52 
0.57 
0.61 



errors are added linearly to the statistical limits. 

The rriEC distribution of the events passing the selection requirements (except for the rriEC cut) 
is shown in Fig. |3| (for NN). For illustration purposes, we have superimposed the Gaussian signal 
contribution and the ARGUS background contribution; both distributions are normalized to the 
results given in Table 0. 



3.4 Maximum Likelihood Analysis 

The likelihood P{f,Xi) for event {i} with a measured set of discriminating variables Xj, and for a 
signal fraction /, is defined as 



P{f,Xi) = fPs{Xi) + (1 - f)Pb{Xi) 



(3) 



The Ps{xi) {pb{xi)) is the product of the normalized signal (background) probability density func- 
tions (PDF) of the individual variables entering the fit: namely, mEC; ^E, m{r]) and the NN or 
Fisher output. The Tyvr invariant mass is not included since the line shape of the oq resonance is not 
well known at present [y]. The data sample is selected within generous sidebands for the variables 
that enter the fit (c./. Sec. |3.1|). The signal distributions are obtained from Monte Carlo simulation, 
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Table 3: Selection efficiencies, expected signal and background yields for the cut-based analysis 
using Fisher or NN. When two errors are given, the first is statistical and the second is systematic. 
We quote upper limits as the main results and compute branching fractions and statistical 
significances, should the observed excess be interpreted as evidence for a signal. The systematic 
errors are discussed in Sec. 0. 





Fisher 


Neural Net 


Signal efficiency (%) 


14.2 


14.6 


Events in signal box 


20 


18 


Events in GSB 


242 


197 


Expected background 
events in signal box 


9.6 ±0.6 ±1.4 


7.9 ±0.6 ±1.1 


Branching fraction (xlO~^) 


8.2l^;[;±1.4 


7.8l|°±1.2 


Statistical significance 


3.0a 


3.10- 


90% CL Upper limit (xlO"'^) 


15.3 


14.5 



> 



20 



15 



Data 

Sig + Bkg shape(overlaid) 

Bkg shape(overlaid) 




5.225 5.250 5,275 5.300 

Mb, GeV/c' 



Figure 3: Distribution of m-EC obtained from the cut-based analysis using the NN. Superimposed 
on the data points are the signal and background rriEC distributions, normalized to the number of 
observed events. 



refined with data from the signal-like charmed B decay mode Dp, and inclusive samples. On-peak 
sideband events, controlled by off-peak data, are used to infer the corresponding background shapes. 
Multi-Gaussian, polynomial functions and cubic splines have been used to empirically approximate 
the reconstructed shapes of the discriminating variables. 

After applying the event selection described in Sec. 3.1, a total of A/" = 9248 candidate events 
without multiple combinations enters the ML fit, corresponding to a signal efficiency of 32.8%. The 
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B(B"^ ao(a,,^ Ti7r)7r) (x 10"^ ) 



Figure 4: Negative log-likelihood function versus the signal branching fraction. 



fit is performed by minimizing the sum 



M 



2lnCf = -2j2^nPif,Xi) , 



(4) 



i=l 



over M events, with respect to the signal fraction /. The fit results in Ng = fj\f = 18.1 ly'4 
signal events for NN, and Ng = 16.2 _y'2 signal events for Fishery. In the following, we restrict the 
discussion to the results obtained with the NN, since it provides the best signal efficiency. This 
choice has been made prior to uncovering the results. The negative log-likelihood function versus 
the signal branching fraction is depicted in Fig. |^. A "toy" Monte Carlo simulation indicates a 
goodness-of-fit of 50%. If the observed signal excess is interpreted as evidence for a signal, the 
corresponding branching fraction is 



BiB"^ 



4(4 



?77r )7r^ 



3.2li;0±l.l) 



X 10 



-6 



(5) 



where the first error quoted is statistical and the second is systematic. The statistical signif- 
icance of this result corresponds to 3.7 standard deviations. The latter is computed using a 
zero-signal toy Monte Carlo simulation. Systematic effects are discussed in Sec. |^. Assuming 
no evidence for a signal, the 90% CL upper limit on the branching fraction, obtained from the 
integral J^ Cf{B')dB' / j;^ Cf{B')dB' = 0.9, is 



B{B^ 



.±/'^± 



7]TT- 



'-)tt^) < 11.5 X 10 



-6 



(6) 



^ From "toy" Monte Carlo studies, the difference of 1.9 events observed between the signal yields with NN and 
Fisher is consistent with the statistical overlap between the samples. Taking into account the 84% correlation between 
the signal yields of the fits with NN and Fisher, the probability to find a difference larger than that observed is 45%. 
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Figure 5: Left-hand plot: distributions of the x variable (R) for data (points with error bars), signal 
Monte Carlo (cross-hatched area) and the background expectation from on-peak sidebands (shaded 
area). The signal region has been expanded in the lower plot. The simulated signal contribution 
is normalized to the ML fit result. A signal excess is observed. Right-hand plot: projection on 
the r]TT invariant mass axis after requiring x > 0.98, removing 44% of the signal and 99.9% of the 
background. 



where the systematic uncertainty has been added linearly. A 65% correlation between the ML fit 
result (^) and the cut-based result (see Table y) is found with a toy Monte Carlo simulation. The 
probability for finding a larger difference than that observed (1.6 x 10~^) is found to be 54%. 
The signal yield of the likelihood function can be displayed with the x-variable, defined by 



Ps -Pb 
Ps +Pb 



(7) 



Its distribution for the data sample used in the ML fit is shown by the points with error bars in the 
left-hand plots of Fig. 0. Also shown are the expectations from signal Monte Carlo (cross-hatched 
area), normalized to the signal yield of the ML fit, and background from on-peak sidebands (shaded 
area). A signal excess and good agreement between data and simulation are observed. 

Since the information from the shape of the rjTT invariant mass distribution is not exploited in 
the likelihood function, one can use projections in order to verify that the signal candidates are 
consistent with the oq hypothesis. We select events with x > 0.98, keeping 56% of the signal while 
retaining only 0.1% of the background. The distributions of the corresponding data events, as well 
as the signal and background expectations are shown in Fig. g. The data are consistent with the 
expected enhancement at the oq invariant mass. 
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4 Systematic Uncertainties 

The main sources of systematic uncertainty originate from the accuracy of the simulation for the 
reconstruction of neutrals, from the tracking efficiency and from particle identification. Dedicated 
studies provide estimates of these effects. A 5% uncertainty is assigned to the r] reconstruction 
efficiency, which represents the dominant error on the selection efficiency. The tracking efficiency 
difference between data and Monte Carlo simulation amounts to 3.9%. An uncertainty of 1.8% is 
assigned due to particle identification. 

Other systematics may arise from the imperfect simulation of the distributions of the discrimi- 
nating variables. They have been studied with control samples. The widths and central values of 
the AE and niEC signal distributions are calibrated with B^ —^ D^{—^ K^tt^)p~ decays, where 
the estimate of systematics using final states with vr'^'s is expected to be conservative since the vr'^'s 
have more combinatorial background than ry's. The r] and ao invariant mass distributions have been 
studied with inclusive data samples. Since the line shape of the ao resonance is not well known Q, 
we have used a generous mass window (c./. Sec. ^^ so that the resulting systematic effect from the 
r/vr mass requirement is small. The Fisher and NN distributions have been checked by comparing 
on-peak sideband data and qq Monte Carlo, where all differences are assigned as systematic errors. 
The uncertainty on the signal efficiency related to the limited Monte Carlo statistics is negligible. 

The contamination from continuum background in the cut-based analysis is evaluated by ex- 
trapolating the number of events found in the GSB, using the known shapes of the AE and ttzec 
background distributions. The consistency of the parameterizations has been checked for different 
regions of the sidebands, and between data sidebands and qq Monte Carlo samples as well as off- 
peak data. Good agreement is found. Associated uncertainties are estimated by varying the shape 
parameters within their statistical accuracies. The same procedure has been applied to estimate 
the systematic uncertainties from the cuts on the event shape variables and the MVA discriminants. 

The estimate of the (mostly charmless) S-background has been performed with generic charm- 
less Monte Carlo. We expect a contamination of 0.3 events from non-resonant vr+vr^vr'^ decays. The 
feedthrough from the unknown decay mode B^ — > aQK~ is determined from the comparison of the 
inclusive a^h' signal yield (i.e., obtained without PID requirements) with the exclusive results of 
the cut-based analysis. No indication for a possible contamination from aQK~ events is found. A 
conservative 5% systematic is assigned to the kaon fraction in the a^i:' signal. 

Table ^ summarizes the relative systematic errors used for the cut-based analysis. Where they 
differ, uncertainties are given separately for NN and Fisher. 

Systematic effects specific to the ML analysis arise from differences in the background shapes of 
the MVA variables between data and Monte Carlo simulation (15%). The systematic uncertainty 
from the correlations between the discriminating variables entering the fit, in particular those 
between the signal distributions of AE and w-eCj amounts to 2%. The total systematic error for 
the ML fit is 18%. 

5 Conclusions 

The preliminary analysis reported in this paper shows an excess over expected background of 
B^ — > aQ'(980)7r^ events, excluding the zero-signal hypothesis at the level of 3.7 standard deviations. 
Interpreted as evidence for a signal, the excess would result in the branching fraction B{B^ — > 
ao'(ag' — > 777r"'")7r~) = (6.2l2'.5il-l) >< 10~^, where the first error quoted is statistical and the second 
is systematic. This corresponds to an upper limit of B{B^ -^ '^t^'^t ~^ r/7r"'")7r~) < 11.5 x 10~^ 
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Table 4: Summary of systematic uncertainties for the cut-based analysis. 



Source 




Uncertainty 


f] finding 




5.0% 


f] mass and resolution 




0.7% 


ao mass and width/resolution 




0.1% 


Track finding 




3.9% 


PID 




2.6% 


AE 




1.7% 


ruEC 




0.1% 


9t and ^tp 




1.0% 


MVAs 


NN 


: 1.0%, FI: 2.0% 


GSB qq background estimate 


NN: 1.3 events, FI: 1.5 events 


Charmless B background 




0.3 events 


OqK' feedthrough 




5% 


MC sample statistics 




0.2% 


B counting 




1.5% 


Total error on branching fraction 


NN: 


15.7%, FI: 17.0% 



at 90% CL. We emphasize the use of improved, linear and non-linear multivariate background 
suppression techniques and optimal signal extraction criteria for rare decay searches. 
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Appendix: Cut Optimization 

The cut on the MVA variable output Xout (the final cut being denoted Xo^^) is applied following a 
criterion designed for rare decay searches. In particular, it does not require one to know (or guess) 
the branching fraction of the signal one is looking for. When searching for a rare signal, one wants 
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to rule out the null hypothesisFI: "there is no signal in the data sample". For this purpose, one 
defines the confidence level 

oo 

CL{N,xZt)= J2VniN,), (8) 

n=N 

where N is the number of events retained by the final cut, A^^ is the expected number of background 
events passing the final cut, and Vni^b) is the corresponding Poisson probability distribution. If 
CL is below a certain threshold {e.g., 5%, meaning that the hypothesis is excluded at 95% CL), 
one excludes the absence of rare decays, to this level. The idea is to adjust a;™^, hence varying 
both signal and background efficiencies to reach the lowest confidence level, on average. 

Cut Optimization for Known Branching Fraction B. First, if the branching fraction is 
known, then, for a given x^^*, one can predict the expected number of Ng signal events. The 
expected value of the CL of Eq. ^ (i.e., the average over a large number of hypothetical experiments) 
is 

oo 

(CL)(xrt) = E T'NiNs + N,) CL(iV, x^) . (9) 

N=0 

The optimal cut on Xout, denoted x°^^{B), is the one which minimizes the above average, i.e., which 
leads to the clearest rejection of the (wrong) hypothesis: "there is no signal in the data sample" 



{CL){xZm = (CL)^i,;,c..(xrt) • (10) 



Figure ^ (left hand plot) shows the confidence level achievable for a given number of expected 
signal events. Also shown are the variations of the optimal cut and the cut obtained when using 
the Gaussian criterion S'^/{S + B) for optimization. 

Cut Optimization for Unkno^vn Branching Fraction. When no expectation for the signal 
branching fraction is available, Eq. ^ cannot be used to define the optimal cut value of x^^^. The 
idea is then to replace the target branching fraction by a target confidence level. One chooses a 
particular value for the target CL (denoted CL*^"*) one is aiming at to define the optimal x'^^ value 
as the value x°^^ (now without {B) as an argument) for which the equality 

(CL)(x:p*) = CL-*, (11) 

is reached for the smallest B. In the following we use CL™* = 0.05. The performance for unknown 
branching fraction is illustrated for the NN discriminant in the right hand plot of Fig. ^. 
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